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ABSTRACT 


Th« feaslbllicy of using range/poinclng angle data such as might be 
obtained by a laser rangefinder for the purpose of terrain evaluation in 
the 10-40 meter range on which to base the guidance of an autonomous rover 
has been further investigated. The decision procedure of the Rapid Estima- 
tion Scheme for the detection of discrete obstacles has been modified to 
reinforce ihs detection ability. With the introduction of the logarithmic 
scanning scneme and obstacle identification scheme, previously developed 
algorithms are combined to demonstrate the overall performance of the inte- 
grated route designation system using laser rangefinder. In an attempt to 
cover a greater range, 30m to 100m, the problem estimating gradients in the 
presence of positioning angle noise at middle range is investigated. 
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DATA ACQUISITION AND PATH SELECTION DECISION MAKING 
FOR AN AUTONOMOUS ROVING VEHICLE 


I. INTRODUCTION 

Current national goaln now include an autonomous rover mission to Mars 
in 1986 and a possible sample return mission for Mars in 1990. The values to 
be derived from these missions will depend on the extent of the ranges of ex- 
ploration which are permittid by the roving devices and the effectiveness of 
their guidance systems. 

The large round-trip communications delay which ranges from about nine 
to forty-five minutes, depending on the positions of Mars and Earth, precludes 
direct earth control except for emergency situations. Systems capable of 
sensing terrain, detecting hazards and selecting safe but efficient paths to- 
wards desired locations are essential. 

In seeking alternative concepts for the unmanned exploration of Mars, 
consideration should be given to: (a) the character of tha rover, and (b) the 

path selection system. 

The mobility, maneuverability and stability of the rover will have a 
direct influence un the availability of safe paths to the desired locations. 

A rover of limited ability to handle in-path or cross-path slopes or boulders, 
craters, trenches, etc. , will be forced to follow tortuous paths to the speci- 
fied destinations. Indeed if its mobility is sufficiently limited, there may 
in many cases be no acceptable path. Tlius, one objective of research and de- 
velopment programs should be the conception and evaluation of rovers of ex- 
ceptional mobility. 

The quality of the path selection systems which are available will also 
have a profound influence on the mission. Tlie term "path selection system" 
is intended to encompass the devices employed to gather terrain information, 
the procedures by which these data are processed and the algorithms used to 
select the optimal path. A path selection system possessing a poor capability 
to discriminate between safe and hazardous features will have to be biased 
conservatively to protect the rover. In such a case, an unnecessarily longer 
path will have to be followed. Indeed, in certain situations, the conservative 
bias may exclude all possible paths and segments of the mission will go unful- 
filled. On the ocher hand, .1 path selection system of high accuracy and per- 
ceptiveness will make it possible to take full advantage of the mobility of 
Che rover. Shorter paths will be followed, less time will be consumed and more 
areas will be explored. 

Research programs at Rensselaer have been aimed at both of these major 
objectives. The research program described herein has been focused on the 
feasibility of using ronge/poinClng angle data, such as might be acquired by a 
laser rangefinder, to characterize the terrain in such detail as required for 
safe designation of route for an autonomous rover and on the development of 
path selection algorithms employing such information. 
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II. SUMMARY OF PROGRESS 


The research program during the past year addressed Che following 

Copies: 


1. Mlddxe range gradlenC esdmaclon wlch noisy eleva- 
vaclon angle. 

2. ObsCacle dececClon by rapid esClmaClon scheme vlch 
modified decision Cree. 

3. InCegraCed paCh selecClon syscem uClllzlng obsCacle 
dececclon and Cerraln modeling algorlchms. 

FlrsC, an approach Co Che problem of esclmaClng In-paCh and cross-paCh 
slopes In Che presence of poslclonlng angle noises aC middle range, 30m Co 
100m, has been developed. The change In slopes has been expressed as a func- 
Cion of Che noise percurbaClons affecClng Che elevaclon angle measuremencs. 

Second, decision procedure of Che Rapid EsClmaClon Scheme and Che evalu- 
aClon of Che probablllcles are modified Co reinforce Che deCecClon ablllCy. 
SlmulaClon resulca show ChaC Che modified version succeeded In reducing almosc 
all Che false deCecClons aC far dlscances. Also, Ic showed Increased deCecClon 
dblllcy ac near dlaCances. 

Third, wlch Che IncroducCion of logarichmic scanning scheme, and obsCa- 
cle Idenclf IcaClon scheme, an overall performance of Che IncegraCed paCh selec- 
Clon syscem is demonscraCed. This Incegral pach selecClon syscem defines Che 
coordlnaclons among Che previously developed algorlchms such as Rapid EsClmaClon 
Scheme, cerrain modeling, and paCh selecClon algorlttims. The compucer simula- 
Clon resulcs are promising In ChaC Che rover Is capable of elecclng a locally 
opclmum path for low risk and energy expendlcure while malnCalnlng a desired 
cargec dlrecClon. 

III. DETAILED PROGRESS REVIEW 

TASK A. GradlenC EsClmaClon aC Middle Range - Ricardo Mediavilla 
FaculCy Advisor: Prof. C. N. Shen 

1) Genesis of Che problem - SCochascically Spaced Rows of a Measuremenc Macrlx. 

The middle range, some 30 Co 100 mecers dlsCance, poses a problem for 
Mars rover navlgaClon. If a laser rangefinder Is locaCed aC Che Cop of Che 
masC of Che rover, say Cwo meCers above Che ground, a horizontal terrain point 
100 mecers away has a very shallow elevaclon angle of about 0.02 radians. A 

point 2 mecers inward from the first point (98 meters from the laser mast) will 

have an elevaclon angle B of 0.0204 radians, which makes an included angle AB 
of 0.0004 radians, or about 1.5 arc-minutes. This approaches the standard de- 
viation of the measuremenc of the elevation angle 5B, which at present can be 

assumed to be 1 arc-minute. In this case the rows of the range matrix displace 

stochastically. 

2) Introduction. 

Terrain iso-gradlents can be used Co evaluate passable Cerraln regions 
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for an unmanned roving vahlcln to h« U5«d on Mars' surface. Range data, 
measured in a spherical coordinate system, is stored in a range matrix. 

Such range matrix consist of a Listing of the distance measurements from a 
ranging device to the terrain surface as a function of asimuth and eleva'- 
tion angle measurements. So far, the range data matrix tias been assumed to 
he fixed xn spacing since the noise in the elevation and azimuth directions 
is small. However, for the shallow elevation angles near the skyline the 
magnitude of the noise in the elevation angle can approach that of the sig- 
nal, displacing the information in the range data matrix. 

3) Scanning scheme 


Range data are stored in the range matrix 
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where 6, the elevation angle, varies with index i, and d. the azimuth angle, 
varies with index j. 


In order to change the data spacing as a linear function of the radial 
distance p. from the rover it is desired that Ap ■ P.A0. (See Figs. 1-a and 
1-b.) ^ 


APj^ - pj^Ad 


(1) 


tand 


■>1+1 


tand 


i+1 p^-Pj^Ad 


p^tandi^l - 


j HAd 

PiAdtandi^^ - r-^ 

Thus, a recursive equation for the calculation of Is given by 

. ► o HA9 

Ap^tand^^j_ - 


( 2 ) 
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Fig. 1-b. Range Measurement. 
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4) CompuCaclon of row spacing. 

Taking Ad as constonc In aquation (2), the alavacion angle is re- 
cursively calculated. Table I shows some typical calculations of the eleva- 
tion angle where equation (2) and a value of Ad ■ .03 radians has been 
used . 


Note that as the distance from the rover becomes smaller, A6>>! arc- 
minute, and the error 36 in elevation angle measurement becomes negligible 
compared to AS, the angle increment. 

Table I 


Recursive Calculation of Elevation 
Angle 6 for Horizontal Terrain 



Distance from 



A6 

Elevation Angle 6 

Mars Rover p 

Ap 

A3 

(arc- 

(radians) 

(meters) 

(meters) 

(radians) 

minutes) 

.02499 

92.11 

1.53 

.00078 

2.69 

.02577 

90.58 




.02910 

83.28 

2.48 

.00090 

3.11 

.03000 

80.80 




.05687 

30.98 

.64 

.00175 

6.03 

.05862 

30.34 


. 

\ 


.15485 

13.45 

.40 

.00475 

16.32 

.15960 

13.05 




i. , , 


S) Error in slope due to error in elevation angle measurement. 

Fig. 2 shows that at far dist.inces from the rover the error in the 
elevation angle measurement can seriously offset the slope calculation. 


Let 36 represent the error in the elevation angle measurement. The 
following quantities are assumed to be: 


0^ ■ 100 meters H ■ 2 meters 
56 ■ 1 arc-minute 


Then from Fig. 2 we have 




and 


’^i+l 


H 

sin(6j^-*-36) 


98.586 


H 


sin(6^^j^-36) 


Noise Effect in the Calculation 


of In-Path Slope at 100 Meters 



3. Noise Effect in the 


Calculation of In-Path Slope at 30 Meters. 
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where hee been calculeCed using equsclon (2). The terrain slope la given 

by 


slope 


l^lSlndi-^1 r^slnJ^ 




(3) 


Fig. 3 shows Chat at a distance of 30 meters from the rover the error In 
the elevation angle measurement does not have as serious an effect as In the 
case of 100 meters distance. Let 


■ 30 meters H ■ 2 meters 
66 ■ 1 arc-minute 

'1 • ; r ac"p;6;) ■ 

''1+1 ' .ln(8j^^-«6) ■ 

Using equation (3), It Is obta5u'.d 

slope ■ .026 


In a flat terrain condition on error of 1 arc-minute In the elevation 
angle measurement will produce an In-path slope of .323 when the distance from 
the rover Is 100 meters, and an In-path slope of only .026 when the distance 
from the rover Is 30 meters. 


6) Evaluation of the In-path slopes from range slopes. 


So far the approach that has been followed Is that of recursively esti- 
mating the range slopes In spherical coordinates, and then transforming to a 
cylindrical coordinate system to obtain the terrain slopes. \a seen from 
Fig. 4 the equations describing the coordinate transformation are given by 


P 


ij 


rijCosBij 


(4) 


9 


iJ 



(5) 


Z 


ij 


H - r^jSlnS^j 


Taking differentials of equations (4) and (6) one obtains 


dZ 


3Z 

96 


d6 + dr 

jr 


dp 


9o 

96 


dS + dr 
9r 


( 6 ) 

(7) 


( 8 ) 
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Th« ln-p«th carrain slopa bacomaa 


:z. IZ . dr 

« , Jx/dg . 38 3r d6 

dp * dp/dS Ifi. ^ l£ . ^ 

38 3r d8 

whara dr/dB la Cha in-path ranga aiopa. 

7) Evaluation of tha In-path alopaa la tarma of tha Invar aa of tha ranga 
alopaa. 

At nlddla ranga tha arror In tha alavatlon angla maaauramant 88 cannot 
ba naglactad. (Saa Fig. S.) 

A A 

Lat r. and r .. ba tha baat availabla aatlmataa of tha ranga naaaura- 
manta. Daflna ^ 



• “i - “i+i ' “i+i - *i 

(10) 


“ • ■ h ■ “ui 

(11) 



(12) 

As 

equations 

before, the spherical to cylindrical coordinate 
are '•■'ven by: 

transformation 


A 

p ■ rcoaS 

C3'» 


e - 8 

(14) 


Z - H - r 

(13) 


Taking dlf ferentlala of aquations (13) and (15) one obtains 

dZ - If dS + •— dr (16) 

3r 

dp • dS + •^ dr (17) 

3r 

Assuming that tha arror In tha range measurement r has a standard 
deviation of 5 centimeters. Its effect can be neglected as compared to the 
effect of error In the elevation. angle measurement 3. So, this time deriva- 
tives are taken with respect to r. 


dp 


dZ/dr 
H /dr 


3Z^>4 

38 dr 3r 

3p 36 2s. 

3 ^ 3 ; 


(18) 


Comparing equation (18) w'.th equation (9) one can see that in equa- 
tion (18) r is considered known while 3 is stochastic. In equation (9) the 
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Fig. 5. Displacements in the Range Measurements Caused by 
Noise in the Elevation Angle Measurement. 
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error in 3 la negligible, and eo d ia conaidered Co be fixed. For ch> ae 
reaaona dif ferenciala are taken vrich reapecc Co 6 in equation (9) onu with 
reapect Co r in equation (18). 

From Che apherical to cylindrical tranafornation equationa, (13) 
through (IS), one obtains: 



3Z 

38 " 

-rcos8 

-*T • cos8 
3r 


3Z 

* 

Dr* 

-sin8 

It ■ 



d8 



(-rccsS) 


-r - tan8 

dr 

dp 

(-rslnd) 

^ cos8 
dr 

-r(tan8) ^ 1 

dr 


(19) 


(i B 

where — r ia Che Inverse of the range slope, 
dr 

* J U 

Le»‘ Con y - -r ^ . Tlien 
dr 


dZ tany - tanp / ,v 

37 ■ 1 t tan>t.na • 


( 20 ) 


When Ar 0 in equation (20) the terrain will be interpreted as a 
vertical wall and Cany ■* creating a problem in Che evaluation of dZ/dp. 

8) Error analysis. 

Fig. 5 dhows how noise 68 in the ^elevation angle measurement 8 can 
displace the range measurement r. Let r represent Che actual range measure-- 
ment after accounting for Che noise in the 8 measurement. Let 68. be 
Oaussian random noise of zero mean and 1 arc-minute standard deviation. <\n 
expression will be derived relating a change in slope due Co 68. 



^1+1 " ‘■i+l®^"‘^i+l 
’*1+1 " *^l+l‘^°“^i+l 


sin(8^^^ - 68^^^) 


sin(8^ + 68j^) 
Hsindj^^j^ 

'*‘"('* 1+1 - 
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% KsinS^ 

>\/ UcoaB. 

X. - r.coaS. ■ . . , — — Va v 

lit aln(3^ + 60^) 


D«fin« 3S as tha change In slope due Co a perturbation In 0^ and 


5S - - ^ 

*1+1 ’ *1 

ain0j^_^j^8in(0j^ + 68^^) - 8in6^sin(0^_j_j^ - 
coa0j^^j8ln(0~ + 50^) - cos0^aln(0j^^j^ - 


( 21 ) 


eln0^_^j^(ain0^co860^ + co80^8in60^] - 8in0^(sin0^_^j^co860^^j^ - co80j^_^j^8in00^^j^l 

co80^^j^(8ln6^co860^ + co80j^8ln(i0^] - co80^[8in0j^^j^co860j^^j^ - co80j^^j^8in56j^^j^] 


Since 58 Is small, Che following approximations can be used: 


( 22 ) 


cos58^ • 

1 

(23) 


1 

(24) 

slndp^ ■ 

501 

(25) 


'«1+1 

(26) 


Approximations (23) through (26) are valid since 58^ Is Gaussian with zero 
mean and 1 arc-minute stcodard deviation. 

ain0^^j^[sin8^ + 56^cos8^] - sin8^[sln8^^^ - ^S^^j^cosB^^j^] 
cos8^_^j^[sln6^ + 58j^cos8j^l - cos8j^[sinBj^^j^ - 


5S c - 


56^tan8^^j^ + 58j^_^j^tan6^ 
tanS^ - tan8j^_j,j^ + 50^^ + 


(27) 


If It Is also assumed Chat 58 > equation (26) becomes 


5s 


56^(canBj^ + canS^^^j^ 

(tanS^^j^ - tan8^) - 260^ 


(28) 


From equaC.^on (10) It Is obtained that 


can8 


1+1 


tan8^ + tanAB 
1 - tanB^tanAB 


tan(B^ + AB) 


(29) 
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If both •mall, Chen 

1 - canS^can/^6 -*• 1 and 


equation (29) becomes 


6S 


6S 


tan6^^j^2i tanfi^ + tanZiS 

56 , 


« (tan3^ + *^“^ 1 + 1 ^ (CanAB) - 266, 


tan6 + tanB, , , 2 v,_ 

- ( — 

1 


Some method of analysis should be developed for obtaining Che variance of 
6S If 66^ Is Gaussian. 

TASK B . Obstacle Detection by Stabilized Rapid Estimation Scheme with Modi- 
fied Detection Tree - C. S. Kim 
Faculty Advisor: Prof. C. N. Shen 

Rapid Estimation Scheme, Reference 1, was developed earlier to detect 
Che outlines of discrete obstacles by processing noisy contaminated range 
measurement data obtained from laser rangefinder. It was found Chat Che sys- 
tem equations of Che previous R.E.S. are unstable, since Che eigenvalues of 
Che system matrix lie outside of unit circle in the 8-plane. For the reasoi: 
above, Che previous R.E.S. generated many false detections In Identifying Che 
outlines of the obstacles at Che ranges greater Chan 30m or less than 10m. A 
stable system model was formulated by Introducing weighting parameters and re- 
versing Che direction of Che recursive estimation and detection procedure. 

As a result, increasing accuracy In state estimates reduced the number of 
false detections. It was reported. Reference 9, that Che stabilized R.E.S. 
succeeded in detecting the complete outlines of obstacles at 40m away from the 
rover. However, It was also pointed out that Che conditional probabilities 
following Che Input estimation were evaluated by approximation. Thus, the 
effect of Che approximation was partially compensated by assigning excessive 
penalties to false detections In Che Bayes decision process. For the reason 
above, the detection ability of Che scheme Is weakened and Che R.S.E. was not 
able Co detect an Input whose signal to noise ratio Is less than 4.5. 

In an effort to Increase Che detection ability of the R.E.S., modlflca*' 
Cions have been made to Che structure of the decision process and the input 
estimation process. Simulation results show Chat this modified R.E.S. generates 
only few false detections In Identifying the complete outlines of obstacles at 
40m without assigning excessive penalties to false detections. Also, the modi- 
fied R.E.S. detected half of the bottom edges of the boulder at 4m where Che 
signal to noise ratio is In the neighborhood of 1.5. The following describes 
modifications made to the decision process and input estimation procedure: 

R.E.S. estimates the state of the system by using an input estimation 
scheme and it detects large inputs of unknown magnitude by utilizing a decision 
tree and the Bayes decision rule. The impulsive input in the system equation 
corresponds to the existence of an edge of an obstacle. Fig. 6 shows modified 
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Fig. 6. Basic Decision Tree with Input Estimations. 




m-Aur; 



Fig. 7. Consecutive Decision Tree. 
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decision crca, each of which reprasanca a particular hypochaals as follows: 

H^; Impulslva input exists between stages i and i-M. 

H21 iapulsive input exists between stages 1-fl and i-t-2. 

H^; inpulsive input exists between stages i-f2 and i-f 3 . 

H.; no input exists between stages i and i- 4 ^. 

U^th reference to Fig. 6, along branch a, the conditional state estimates x^^ 2 > 
x^.f3 and x^>4 are obtained by using Kalman filter equations, since assumes 

no impulsive input between stages i and l*f 4 . The procedure to estimate the state 

at stage i-fl, will be discussed later. 

For branch 3 , the state estimates for stages i-fl and i-t -2 ar^ the same as 
branch a. Since H3 assumes in input between stages i-f 2 and 1 *^ 3 , X14.3 is obtained 
by utilizing an input estimation scheme. As shown in Reference 2 , the input es- 
timation 0(4.2, obtained from 2^4.4, since the existence of an input 014-2 does 
not affect Z14.3 but Zi4^ directly. As a result, X14.3 is obtained by utilizing 
measurement data z^w one stage ahead. Since Z14.4 has already been used for the 
estimation of X14.3, it cannot be used twice. Thus, the following stage estimate 
X14.4 is obtained by prediction only. As is noticed from above discussion, we 
need at least 4 stages to distinguish branch 3 from branch a. 

For branch 2 , the s^ate estimate at stage i>l is the same as X]^4-i. How- 
ever, the stage estimate xi4-2 is obtained by using an input estimation scheme, 
p^ls is because H2 assumes an input between stages i-t -1 and i-f 2 . As a consequence, 
^ 1+2 is^estlmated by using measurement data Z14.3 one stage ahead. The following 
itate X14.3 is estimated by prediction only, since Z14.3 has already been used for 
che estimation o^ xi4.2« To make branch 2 comparable with branch a and branch 3 , 
stage estimate X14.4 is obtained by using an input estimation scheme, since 
assumes an input between stages i and i+ 1 . The following state X14.2 is estimated 
by prediction only, since Z14.2 has already been used for the estimation of x|^4.i. 

To make br^ch 1 to be comparable with the rest of the branches, state estimates 
x( 4.3 and x(4.4 are obtained by Kalman filter equations which use the measurement 
data z^ 4>3 and z^4^ respectively. 

As shown in Fig. 7 , which will be discussed later, the state estimate x^4.]^ 
is obtained by either prediction or Kalman filtering. If the state has been 
obtained by utilizing input estimation scheme in the former decision tree, x(4.i 
is obtained by prediction only, since the data Zj^4.j^ has already been used in the 
former stage. On the^other hand, if x^ is obtained by Kalman filtering in the 
former decision tree x(4.j^ is obtained by using Kalman filtering. The conditional 
state estimates along the decision tree in Fig. 6 and the corresponding covari- 
ances are used for the computation of the Joint conditional probability of the 
residues for each branch of the decision tree as follows: 


*^ 1+1 

m 

*^ 1+2 


^+2« ^1+3* ^44!^’ V 

Zi_j_l-E{H f( Xj^ for 


*i+2"^^^ * 1 + 2 ' *1+1^^ " ^i+2"^ ^1+1 *1+1 


m ■ 1,2, 3 , a 
for m - 1 , 2 , 3 a 


(a) 

(b) 

(c) 
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^l-t-3 " *1+3'^^” *l-*-3l*l+2’ " *l-*-3“** *i>2 


*^1+3 ^ *1+3*^^“ *1+31*1+2’ “1^ " *1+3"** ^1+2 *1+2 


(•) 


^+4 • *l+4-^^» *l+4l*^l+3’ " *l+4’» ^1+3 *J+3 

*^i+4 ^ *1+4"^^” fof “ • 1.2 (g) 


1-^4' 1+3’ m' 1+4 " ‘1+3 1+3 


The Joint conditional probabilities are used to evaluate the Bayes risk Index 

ned below: 

**m^liii^r^’^l+l’’^l+2’’^l+3’^l+4l*l’”iB) 

y for H . 
m 

; cost associated with accepting H. when H Is true. 

LJtn t n 


for the branch t as defined below: 


R, ■ E 

m«l,2,3,a 


where P ; priori probability for H 

m D 


If R2 « R ^ or R ^ Is minimum. It Is concluded that there Is no Input between 
stages 1 and 1+1. As a consequence, the state estimate Is accepted as a 

correct state estimate and branch 1 Is discarded. As shown In dotted line in 
Fig. 7, with x°^]^ as a new Initial state, branches 2, 3 and a are extended to 
stage 1+5 and become branches 1, 2 and a respectively. Branch 3 Is newly 
generated and the algorithm proceeds with a new set of data. 

If minimum. It Is concluded that there Is an input between stages 

1 and 1+1. accepte(i as a correct estimate for stage 1+1 and totally 

new set of decision tree Is generated to the stage 1+5 with as a new 

Initial state. This procedure Is Illustrated In chain dotted line In Fig. 7. 

It should be noted again that should be obtained by Kalman^f liter equa- 
tions when Is accepted as a correct state estimate. When Is 

accepted as a correct state estimate, x^^2 obtained by prediction only since 
the measurement data already been used for the estimation of xi+i In 

the former decision tree. 

Now, the joint conditional probability In equation (2) Is computed as 
follows : 


P 

r 


. ui 

^*’ 1 + 1 ’*^1+2 


m m 
*^l+3’^l+4 




*l’*l+l’*l+2’^l+3’“m^* 


**..^*■ 14 . 3 1 * 1 ’* 


m 


‘1+2’ m 


,H.) • P,(r 


m 


1+2 


X . 9 z 


1+1 ’“m^ 


^(^+1 * i’V 


It should be pointed out that the last term of equation (12) , 

Pj.(r. -Ixj^.H^) Is the same for all branches, since rj+j^ Is common to all the 
branches as can be seen In the decision tree of Fig. 6. ^For this reason, 
reduced Joint conditional probability Pr(*‘“+2 *l+l.^n will be used 

for comparison among the hypothesis. 

.Vs Markov property has been preserved in the procedure of state estima- 
tions, the reduced Joint conditional probability for each branch is written 
as follows: 
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i 1*1 


1 1*1 


Pr^’^1^.2l*i+l'“l^ 


2 2 2 


**r^^H-2’*^l-»-3’'^t-f4l*i+l’‘*2^ " **r I *l-*-3 '**2^ ’ **r^‘^l-*-3l*l-*-2’*^2^ 


(Aa) 


**r^'^i-*-2l*^l-*-l'**2^ 


^r^'^l+2’'^l+3’^l-H4l*i+l’‘*3^ " **r^*^l+4l*l+3’^3^ * **r ^'^l-*-3^*l+2’‘*2^ 


(4b) 


P (•■** lx* H ) 


(4c) 






-w“.> 








I‘l+l-"a> 


(4d) 

For Che comDuCaClon of Pr (>^14-2 1 1° equaclon (4a), Che expecced 

value of residue ri4.2 is HPj^^j^Bu^ as shown in equaclon (6) of Reference 2. 

As a resulc, Che correcC evaluaCion of Che probablllcy would be: 


where 


and 


^<^i+2l*?+r«l> - T T' i"V """ 

*‘"•''1+2 1 

1 .1 T 1-' .1 

*^1+2 " "i+2’''i+2’’^l+2 

w [^2 - ^^ U >* 

"i+2 * "i+2 ” ^^1+1 ®"i 


2 1+2' 


(5a) 


(5b) 

(5c) 

(5d) 


Since u. In equaclon (5d) Is unknown, Che besc approxlmaclon would be 
Co subsClCuCe ^^i+lB“|’ ^^<4.1®“} which Is derived In equaclon (c8) of 

Reference 2. 

*l 

Thus, Che residue h^<^°<aes: 



1+2 


- HF! 


1+1 


Bu. 


*■1+2 


«^i+l 



(5e) 


Since Che InpuC esclmace uj^ Is obcalned from one ihoC esclmaclon by using 
2^4.2, 1C does noC dlsclngulsh noise from signal. This approxlmaclon will resulc 
In a decermlnlsclc quanclcy as shown In equaclon (5e) , and exacc evaluaclon of 
Che probablllcy Pr ^*’1+2 1 ’'i+l»hl) 1® possible. Ic is^found ChaC similar prob- 
lems30ccu^ wlch Che evaluaCion of Che probablllcles Pr (>^i+3 1 ^*4.2 .H2) , and 
®r (>^i+4 1 *i+3»h3^ 1“ equacions (4b) and (4c) respecCively. Thus, Che j.c.p. for 
each branch becomes: 
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Branch i: ‘ %^*'l+3l*U2'**l^ 


Branch 2: 


% ^ 1+4 * ' 1+2 ' * 1+1 ’ i+3 ’ **2 ^ 


**r^''i-t4l*i+3’**2^ ’ 


i+2l*i+-l*“2^ 


(6b) 


Branch 3: ^p^^£+3»*^i4.2^*l+l*^3^ " ^r^^ 1+3 ^*1+2 ’**3^ * ^r^*^l+2 ^*1+1’^3^ (6c) 


For branch a, all cha conditional probablllclas In aquation (4d) can 
ba avaluatad. Slnca tha saaa number of probabilities are needed to be com- 
pared with other branches, only first two probabilities are chosen to be used. 


>r«>ch a: - P/V3I Vj.",.) 


^r^’^l+2^*l+l’**o^ (W) 


^Cl 

4s^l^ noticed In equations above, the estimation of the states 
and are not needed for the evaluation of the probabilities. 

In the procedure of redefining the J.c.p. of each branch as above, the J.c.p. 
of U3 and H become the same, since the two hypotheses share the same branch 
till Che stage 1+3. As, a result, branch 3 Is deleted from the decision tree 
and the state estimate In branch a Is not needed. The final structure 

of Che decision tree with the residues to be compared are Illustrated In Fig. 8. 

Now we can explain heurlstlcally how the R.E.S. distinguishes a large 
Input Co Che system from a large observation noise. For Che purpose of expla- 
nation, assume Che diagonal element of the cost matrix is zero and all others 
are equal. Also the priori probabilities for Che branches are equal. It can 
be shown that this leads Co a maximum likelihood decision rule for multiple 
hypotheses except that the R.E.S. proceeds one stage for each decision recur- 
sively. Since each of the probability functions is Gaussian, the one with 
smaller residue will have larger probability. One may analyze the effect of 5 
cases on the residuals for the 3 branches. The 5 cases are: 

1. A large impulse input u^ between stages i and i+1. 

2. No input, but large observation noise with 

3. A large impulse input between stages i+1 and i+2. 

4. No input, but large observation noise with 

5. No input, but large observation noise with 

.1 

For case 1, there is an input between stages ^i and i+1; then would-be a 

correct estimate and the following residues rj^.f3 in equation (le) and ri.t4 In 
equation (Ig) becomes small. On the other hand, x^^.. which does not include 
input estimation of uj^, would be a wrong estimate. Large residues result in 
^i+2 equation (Ic) and r"^^^ in equation (Id). For case 2, there is a large 
observation noise with ^i+1 includes a spurious input estimate which 

is not present and Che following residuals r[+3 and r^^^ become large. On the 
other hand, x^.t-1 is obtained by filtering the observation noise by Kalman filter 
and Che following residuals r^^2 *^1+3 will be small. This line of reasoning 

leads to Table II which illustrates Che qualitative picture of Che residuals for 
Che five cases mentioned above. 
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Fig. 9. Edge of Boulder (left) and Crater (right) 
at SOm, outlined by R.E.S. 
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Table II 

caaa 1 case 2 case i case 4 caae 5 


large large large large large ^ 


Branch 1 r^^^ 

small 

large 

large 

large 

smai.x 

'U 

small 

large 

large 

small 

large 

Branch 2 r ®^2 

large 

large 

small 

small 

small 


large 

small 

small 

large 

large 

Branch r ^^2 

large 

large 

small 

small 

small 

a 

^i-»-3 

large 

small 

! 

large 

large 

small 


Looking at Table II, one cai notice that: 

In case 1, branch 1 is Ikely to have the largest probability, 
and detect th presence of Input. 

In case 2, branch 1 will t.xve the smallest probability and the 
scheme will succeed to Ignore the measurement noise 
proceed to the next stage. 

In case 3, branch 2 will havs the largest probability. The 

scheme will proceed to the stage 1-fl and detect the 
Input afterwards. 

In case 4, all the branches are equally likely, since all the 

branches have one large residue and one small residue. 

In case 5, branch a Is most likely, and the scheme will proceed 
to the next stage. 

It has been shown qualitatively that R.E.S. would work very well except 
for case 4. Since all of the 3 branches are equally likely In case 4, a 
slight biasness In the cost assignment or priori probability will handle the 
situation. 

Computer simulations were performed for a boulder-crater obstacle pair 
located 40m and 4m away from the rangefinder. Range data were obtained from a 
previously developed program which simulates the rangefinder measurements. 

Tne measurement data were corrupted by a Gaussian noise with Its standard de- 
viation 0.05m. The rangefinder was assumed to be located ?ir above ground, and 
the boulders and craters were of hemispherical shape with 2m in diameter. For 
the cost matrix, the diagonal elements Cxi, ^22 ^33 assigned the value 

zero, since correct decision should cost the least. All the off diagonal ele- 
ments are set to 1.0, so that no excessive penalties are assigned to the false 
detections. In the simulation with obstacles at 40m, elevation angle 3 of the 
rangefinder changes from 0.7° to 3.3° with constant increment A8 ■ 0.04° and 
azimuth angle changes ftom -5° to 5° with A ■ 0.2°. The simulation result 
in Fig. 9 shows that the modified R.E.S. succeeded in reducing almost all the 
false detections without using excessive penalties to false detections. In 
the simulation with obstacles at 4m, S changes from 8.21o to 45.52° with 





Ad ■ 1.49^ And Ch« axlauch angle 6 changes froa -45° Co 45° with ■ 3.6°. 

Fig. 10 shows Che ouCline of obscacles excrscced by aodifled R.E.S. Ic Is 
noced chac Che sodlfied R.E.S. dececced half of Che boccoa edges of Che 
Doulder ac 4« tfhlch could noc be dececced by Che previous R.E.S. Approxl- 
mace signal Co noise raclos aC Cha edges of Che obscacles were coapuced ac- 
cording Co Che foraula SNR • u/o , where u Is Che aagnlCude of Che lapulrlve 
InpuC occur Ing ac Che edge of oblcacle and o, !■ che scandard devlaclon of 
Che noise. The SNR aC Che boccoa of a boulder ac 4a la found Co be In che 

neighborhood of 1.4 while che SNR ac Cha boccoa of a boulder aC 40a la abouC 

7.0. This Is che reason why Che dececclon of che boccoa edge of Che bo<ilder 
ac 4a has bean a harder cask Co be perforaed. 

In suaaury, for Che parclcular sec of scace equaclons and aeasureaenc 
equaCloQS chaC we considered, che exlscence of an Inpuc does noC efface che 
nexc aeasureaenc buC effaces Che aeasuieaenc decs afcer che nexc. For Chla 
reason, cha flrsc secclon of che decision Cree becoaas coaaon co every branch, 
chua Ic fs excluded froa coaparlng che j.c.p. The aagnlCude of an InpuC Is 

esclaaced by second residue which Is one shoe esclaaclon. This one shoe es- 

Cimsclon cannoC dlsClngulsh che signal froa noise. The condlclonal proba- 
blllcy for Che residues which Include che InpuC esclaaclon becoaas deceralnls- 
Clc and chesa probablllcles are deleced froa coaparlng che j.c.p. for Che 
branches. To Case whecher che R.E.S. can distinguish a large IcpuC froa a 
large noise, che residues of che branches are coapared quallcaclvely as a 
aeasure In deceralnlng where cha decision scheae Is correcC. 

In conclusion, Che scablllzed R.E.S. wlch modified decision Cree 
generaCed only few false deCecClons In IdenClfylng Che coaplece ouCllnes of 
obscacles wlchouC using excesslvu penalcles co che false deCecClons. Alsu, 
wlch che reinforced deCecClon abllicy, R.E.S. deCecCed half of Che boccoa 
edges of che boulder aC 4m. 

Thus far, R.E.S. assumed chaC each column or row of Che range macrlx 
Is Independenc and Ignored Che InCerdependency between each oCher. Ic is 
thought chat Che dececclon abllicy of che R.E.S. will be Increased If R.E.S. 

Is expanded Co IncorporaCe che InCerdependency between the columns and rows. 

TASK C. Integrated Path Selection System for Martian Rover Usi n g Laser 
Rangefinder - C. S. Kla 
Faculty Advisor: Prof. C. N. Shen 

There have been separate studies on obstacle detection. References 1, 2, 
terrain modeling and evaluation. References 3, 4, and path selection algorlchas. 
References 4, 6. With Che Introduction of Che logarithmic scanning scheme and 
obstacle Idenclf IcaClon scheme, an overall performance of che Integrated path 
selection system for autonomous Martian rover Is demonstrated In this work. 

Wlch reference to Fig. 11, a laser rangefinder Is used co obtain che 
measurement data. A logarithmic scanning scheme Is designed so that Ic can 
provide che obstacle detection and terrain modeling algorithms with a set of 
measurement data meaningful in Che sense of data spacing over Che entire 
scanning area. The design of che scanning scheme becomes complicated due Co 
Che noise elements In Che range as well as In Che positioning angles. By 
processing Che range measurement data obtained. Rapid EsCimaClon Scheme, 

Reference 2, detects Che possible edges of che discrete obstacles. The result 
of R.E.S. Is fed Into an Identification scheme in which che whole set of 
detected points are classified Into subsets. Each subset Is recognized either 
as a boulder, a crater, cr a ridge. Now, the location and sizes of the obsCacles 
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ar« us< ' Co d«C«rmln« Ch« longchs of the paasabl* corridors for Cha path 
aalacclt schama. In parallal with Cha R.E.S., cha InpaCh and cross-path 
slupas at Cha data points ara astiaatad by racursiva slopa asclioatlon schaaa, 
Rafaranca 4. This racursiva algorithm Is not valid at Cha adgas of obsta- 
clas and rldgas whara dlscrata changas In range occur. Since the edges of 
obstacles aid rldgas are known from the obscacla detactlon schama, the slopa 
estimation scheme will Initialize Itself at those points. Inpath and cross- 
path slopes computed at the data points ara utilized to astlmate the terrain 
varlablas at the spina and track points along each corridor. Tha terrain 
variables considered ara Inpath slopa, tilt slopa, wheal deviation and body 
clearance. Since the spina or track points do not coincide with cha measura- 
msnt data points, tha Inpath slopes, cross-path slopes and elevations at tha 
spina and track points ara astlsiatad from two-dimensional Harmlta Interpola- 
tion function. Reference 8, In a rectangular coordinate syst< a. Now, each 
corridor Is evaluated In terms of the terrain variables. The corridor which 
Is expected to Incur the least amount of energy consumption and risk Is 
chosen by dynaoilc programming method as Che optimum corridor In the path 
selection scheaie. 

The range measurement data Is obtained by changing the elevation angle 
and azimuth angle 9j In spherical coordinate system. The three coordinate 
systems used In this work are depleted In Fig. 12: 

spherical coordinate - (1, 0, B) 

cylindrical coordinate - (p, 9, Z) 
rectangular coordinate - (x, y, z) 

It Is desired for the scanning scheme to covur a horizontal terrain 
from 5m to 80m In radial direction from the rangefinder 2m In height. The 
resolution of the rangefinder Is limited by the noise elements in positioning 
angles as well as Che range error. It Is assumed that the standard deviations 
of the elevation angle and the measurement range are arc-mlnuCe and 

a - O.OSm respectively. Scanning ncheme is evaluated°by three criteria as 

r">i •. 


(1) signal to noise ratio of the Incremental elevation 
angle, 

(11) horizontal data spacing In radial direction, Ap^; 

(111) signal to noise ratio of the Incremental range, 

Assuming that B. 's and m. 's have Gaussian distribution, the quantities 
..11 ao 

and a, become: 

Am 

O-o ■ N2 a ■ 1.414 arc-min. 


a, ■ N2 a ■ 0.07m 
A m 


It Is noticed that to have 46^/o^g and Am^^/Og greater than unity, the corres- 
ponding A6. and Am^ should be greater than 17414 arc-minutes and 0.07m 
respectively all over the scanning area. Since the obstacle detection scheme 
needs at least four data points. Inside or on Che obstacle to be Identified, 


Che horizontal data spacing Ap^ determines Che size of obstacles that can be 
Identified. If constant Increment of radial distance Ap > 0.5m Is used, the 
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Che obstecle dececclon icheme can detect obstacles of 2m in size for all 
range. However, AB. becomes smaller at far dlstancss and at o, ■ 80m, 

AB^ becomes 0.54 arc-min. , which is much snuiller Chan its standard devia- 
tim ■ 1.414 arc-min. Even if constant in''rement of radial distance 

Ap^ • Ao ■ 1.82m will result in the incremental elevation angle AB^ ■ 2 
arc-min. at ■ 80m, the obstacle detecclon scheme cannot identify obsta- 
cles »hose dimensions are less than 7m in size. 

To design a scanning scheme which will make Ap^ at near distance as 
small as possible while keeping AB./o^g and greater Chan unity, the 

following variable increment of elevation angle is suggested. With refer- 
ence to Fig. 12, Che relation between Ap^ and AB^ is derived as follows: 

2AB^h/sin2Bj^ ■ Ap^^tanB^^ (1) 

where 0<B^<n/2 and 0<AB^ and h is Che height of the rangefinder. 

If Che quantity Ap^CanB^ in equation (1) is set to a constant, Ap^ 
becomes smaller and oB. becomes larger as B^ gets larger. From this reason- 
ing, Che scanning scheme will result in a Idrge data spacing at far and 
smaller data spacing at near distances, by setting the radial data spacing 
Ap. to equal the horizontal data spacing AB^, it has been derived in Refer- 
ence 7 chat Che quantity Ap becomes smaller logarithmically as the index i 
increases. In Table III, the resulting data spacings and Che corresponding 
incremental elevation .ingles are listed. 


Table III 


Data Spacing and Signal to Noise Ratios 


*^1 

80m 

— 

50m 

25m 

5m 

■'“l ■ ■‘'*1 

2.4m 

1. 5m 

0.75m 

0.15m 


34.3 

1 

21.4 

10.7 

2.1 

ABi 

2.7 arc-min. 

4.3 arc-min. 

8.4 arc-min. 

36.5 arc-min. 


1.9 

3.0 

5.9 

25.8 


From Che measurement data obtained in this manner, R.E.S. can identify the 
obstacles whose dimensions are as small as 0.6m at * 5m and dimension of 
6m at ■ 30m. The range data obtained according to this logarithmic 
scanning scheme are stored in Che form of range matrix 


M - 






( 2 ) 


*'^1 \2 ••• \n 
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To cov*r a^scc^or in horlzoncal remain defined by 5m p - 80m and 
-0.54 rad ■9-0.54 rad. Che required number of rows and columns of Che 
range macrix in equacion (2) are k - 102 and n ■ 37 respecCively. 

Now, R.E.S., Reference 2, processes Che measuremenc daCa obCained 
from logarichmic scanning scheme Co deCecC Che possible edges of obscacles. 

To show Che dececcion abiliCy of • he R. E.S., range measuremenc daCa are 
obCained from a previously developed range measuremenc simulaCion program. 

The general cerrain model is a wavy surface wtch placeaus and valleys ana- 
lyCically defined in recCangular coordinace syscem as: 

z - 0.5 cos(^) sin (^) (3a) 

where Che unics of x, y and z are meCers. Also Che following feaCures are 
super-imposed on Che general Cerrain described by equacion (3a) 

(1) heml-elllpsoldai boulders ac (x, y) ■ (0,17) and (11,34); 

(li) heml-spherlcal boulder ac (x, y) - (1,7); 

(ill) heml-spherical craCers ac (x, y) - (7,45), (-9,44), and (-^*,67). 

The Cop view of Chls Cerrain is shown in Fig. 13. The rangefinder 
is locaCed aC (x, y, z) ■ (0,0,2) and Che scanned area is defined by 
S ■ 0.3805 rad, 8 " 0*025 rad, 9 - 0.54 rad, and 9^ - -0.54 rad. 

rrftie process of geS^aCing Che measuremenc daCa, whiCe Gaussian noise wich 
0.05m sCandard deviaCion is added Co each of Che range measuremenCs. 

The measuremenc daCa obCained are processed by R.E.S. The verclcal 
edges of Che boulders, craCer and ridges are deCecced as and Che side 

edges are deCecCed as 'O' in Fig. 14. In Fig. 15, Cop of Che ridges are 
Idencified by chains of several column deCecClons. Boulders are idenCified 
by consecuClve column deCecCions ac Che cop and single column deCecCions aC 
Che boCCom. CraCers are idencified by single column deCecCions aC far edges 
and consecuClve column deCecCions aC near edges. Ic should be noCed ChaC 
Che use of R.E.S. is subJecC Co produce false deCecCions as well as Co miss 
some of Che deCecCions. 

Ai has been menCioned previously, iC is needed to develop a scheme 
which will IdenClfy Che whole seC of deCecClon polnCs inco subseCs of deCec- 
Cions. By uCillzing Che characCerlsCics of Che boilder, craCer and ridge, as 
well as Che required relaclonshlp becween Che detected points Co form on out- 
line of an obstacle, an algorithm has been developed. Reference 7, to identify 
a sec of detections as a boulder, a crater or a ridge. By using Che idenclfl- 
caClon scheme, Che whole set of detection points obtained from R.E.S. in 
Fig. 14 is characterized either as a boulder, a crater, a ridge cr a false 
detection. As shown in Fig. 15, Che scheme clearly identified each of the 
boulders, craCers and Che ridges. This information is fed into path selection 
scheme in determining Che lengths of the passable corridors. 

Recursive slope estimation scheme, Reference 4, was developed earlier 
to esCimaCe the inpaCh and cross-path slopes at the data points by utilizing 
Che range data. As has been menCloned earlier, this scheme cannot be applied 
where discrete jump in range occurs, such as edges of obscacles and ridges. 
Since Che edges of obstacles and ridges are Identified by R.E.S. and obstacle 
idcnclf icaCion scheme, Che recursive slope estimation scheme is modified to 
reinlclalize itself at those points. For simulation of Che scheme, Che range 
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Fig. 14. Detected Edges by R.E.S. 



Fig. 15 ; Outlines of Obstacles and Ridges 
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daca obcainad for cha Carraln in Fig. 13 ^ra procasaad by tha racurslva 
alopa aaciaatlon schema. To display tha Informaclon, tha gradlant for each 
daca point Is classified to lie within a specified Interval for mapping 
purpose. Each Interval and the corresponding character Is defined In 
Table IV. 


Table IV 

Isograd lent Regions 


Gradient 


Character 

min 

max 

A 

0.0 

0.02 

B 

0.02 

0.04 

C 

0.04 

0.06 

D 

0.06 

0.08 

E 

0.08 

0.1 

F 

0.1 

0.15 

G 

0.15 

0.2 

H 

0.2 

0.3 

I 

0.3 

0.4 

J 

0.4 

00 

U 

ut.de fined 

gradient 


Fig. 16 shows as Isogradlent map describing Che terrain In x -6 plane 
projection. It Is noted Chat there are circular Isogradlent regions on the 
faces of slnusoldel dunes and valleys. Also, high gradients such as 'J' 
are assigned on the faces of boulders and Inside the crater. It Is also 
noted that hidden regions exist at Che back of a boulder or Inside the crater 
where laser beam could not reach. Undefined gradient, represented by charac- 
ter 'u', occurs when Che recursive algorithm cannot be applied due Co a large 
discrete change In ranges between Che adjacent measurement data. The esti- 
mated InpaCh and cross-path terrain slopes and the Iso-gradlent map are used 
for the evaluation of Che passable corridors. 

So far, we have obtained Che outlines of obstacles and slope estimates 
at the data points by utilizing logarithmic scanning scheme, R.E.S., obstacle 
Identification scheme, and recursive slope estimation scheme. Now, with Che 
Informaclon obtained above, path selection scheme is utilized Co determine an 
optimal corridor which minimizes certain performance Indexes. The path selec- 
tion scheme used here Is based on the method outlined In Che previous work. 
References 5, 6. Oue major modification made to the previous work Is the 
evaluation of terrain variables. Since we have already obtained the slopes 
at Che data points, the terrain variables (inpath slope, tilt slope, body 
clearance and wheel deviation) are estimated by employing two-dimensional 
Hermlnte Interpolation function which utilizes Che estimated terrain slopes. 

A detailed description of Che modifications is found in Reference 7. The 
simulation of the path selection procedure was performed for Che terrain illus- 
trated in Fig. 13. With tht. starting location (x, y) ■ (0,0) and the target 
location (0,350), the rover determines the corridor lengths and evaluates the 
terrain variables by utilizing the algorithms described earlier. Finally, 
optimum corridor is chosen by path selection scheme. The cost associated with 
each corridor is tabulated in Fig. 17 where the index (1, j) corresponds Co 
that of the primary corridor CR^ and secondary corridor CR^j . The length of 
primary corridors CR 3 and CR 4 become short due to the presence of obstacles 
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Fig. 17. Corridor Costs.' 
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•c near dlsCancee. The value 100 is assigned to penalize short corridors. 

The teralnal costs for CRj^ and CR3 become larger because of the target de- 
viation and low elevation at the end of the corridor. As Is noted In Fig. 17, 
the scheme selected the corridor of 25m In the direction 9 ■ 0.18 rad for 
the corridors (2,2). In the next path selection procedure, the rover scans 
the terrain at the position (x, y) ~ (-4.5, 24.6). The measured data are pro- 
cessed by R.E.S. to detect the possible edges of obstacles. Fig. 18 shows the 
outlines of craters detected. Crater #1 corresponds to the crater at (x, y) “ 
(-9,44), which was not detected In the first path selection because It was 
hidden by the ridge. Crater #2 corresponds to the crater at (x, y) * (-4,67). 
It was not detected In the first path selection procedure since It Is too small 
compared to the data spacing at the far distance. Now, the rover utilizes 
again the terrain evaluation and path selection scheme described earlier. In 
this second path selection procedure, the corridor of 25m In the direction 
6 ■ 0.18 rad Is selected as optlm<un corridor. The route resulting from the 
two path selection procedure Is Illustrated In dotted line In Fig. 13. 

In conclusion, this work presents an Integrated study of logarithmic 
scanning, obstacle detection, terrain evaluation, and route designation prob- 
lem. The resolution provided by the scanning scheme Is appropriate for the 
detection of small obstacles at near distances as well as large obstacles at 
far. The outline of obstacles and ridges are Identified from the detected 
points. This Information Is used In the design of passable corridors In that 
the corridor lengths become adaptive to the obstacles Identified. The terrain 
variables along each corridor are computed by using the terrain slopes esti- 
mated recursively. The path selection procedure Is made more practical by em- 
ploying a finite dimension vehicle rather than a point mass. The Integrated 
system for the autonomous navigation of Mars rover Is completed by establish- 
ing the coordinations among existing algorithms as well as making modifications 
to them. Computer simulation results are promising In that the rover lu capa- 
ble of selecting a locally optimum path for low risk and energy expenditure 
while maintaining a desired target direction. 

The following are the areas in which more research work is needed for 
Improved performance of the route designation problem for autonomous Mars 
rover ; 

(1) Adaptive scanning scheme. 

It Is needed to develop a scheme which will detect the skyline 
and make the scanning scheme adaptive to the upgraded hills as 
well as downgraded valleys. 

(11) Improvement on the obstacle Identification scheme. 

The obstacle Identification scheme outlined In Reference 7 Is 
Just a beginning of the contour extraction problem. A more 
rigorous study Is needed to employ all the possible Information 
in constructing the outlines of the obstacles and ridges. 

(Ill) Improvement on path selection scheme. 

In the previous path selection scheme. Reference 5, all the 
terrain variables computed are treated In a deterministic 
manner. It Is needed to modify Che previous scheme so Chat 
the stochastic nature of the terrain variables are considered 
in the evaluation of the cost for each corridor. 

Also, in the previous scheme, the relation between the proba- 
bility of impassablllty and the amount of risk Involved Is 
not clear. 
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